VIBRATIONAL MODES OF CIRCULAR FREE PLATES UNDER 

TENSION 



L. M. CHASMAN 

Abstract. The vibrational frequencies of a plate under tension are given by 
the eigenvalues uj of the equation A^u — tAu = uju. This paper determines 
the eigenfunctions and eigenvalues of this bi-Laplace problem on the ball under 
natural (free) boundary conditions. In particular, the fundamental modes — 
the eigenfunctions of the lowest nonzero eigenvalue — are identified and found 
to have simple angular dependence. 



1. Introduction 

The importance of the disk for physical isoperimetric inequahties has been rec- 
ognized since 1877, when Lord Rayleigh [T7| conjectured that the first eigenvalue of 
the Dirichlet Laplacian on a planar domain (the fundamental tone of a fixed mem- 
brane) is bounded below by the first eigenvalue of the disk of the same area. This 
conjecture was later proved by Faber [7] and Krahn [9j[l0]. The disk is also the ex- 
tremal case in the Kornhauser-Stakgold upper bound for the Neumann eigenvalue 
(free membrane), proven by Szego [HI [19] and Weinberger [22] in the 1950's. 

The disk further gives extremal cases for vibrating clamped plates. Plate vi- 
brations are governed by the bi-Laplacian operator just as the Laplacian governs 
vibrations of membranes. The collective work of Szego [18], Talenti [20], and Nadi- 
rashvili [13j [M] (and Ashbaugh, Benguria, and Laugesen [2 |4j [3] in higher dimen- 
sions) established a sharp lower bound for the fundamental tone of a clamped plate 
— that is, the first eigenvalue of the bi-Laplacian with boundary conditions u = 
and du/dn = 0. As in the Faber-Krahn inequality, the disk provides the lower 
bound on the first eigenvalue. 

In this paper, we consider the free (unconstrained) plate problem, with the goal 
of identifying the eigenfunctions and eigenvalues of the ball, and in particular the 
fundamental mode. We consider the more general situation of a free plate under 
tension and factor the eigenvalue equation in order to show the ball eigenfunctions 
can be written in terms of Bessel and modified Bessel functions and spherical har- 
monics. We identify the fundamental tone using the plate Rayleigh quotient and 
the specific form of the free boundary conditions. 

Free boundary conditions are more complicated than clamped ones. This makes 
determination of the free circular plate eigenfunctions somewhat more difficult than 
for the clamped circular plate [11] Chapter 2], which has only been treated in the 
zero tension case as far as we know (see, e.g., the survey paper [3]). Other boundary 
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conditions exist, such as the hinged plate investigated by Nazarov and Sweers [15j . 
and the simply supported plate. 

The results of this paper identifying the fundamental mode of the unit ball are 
essential to the solution of the free plate isoperimetric problem, which we present 
in [5] . That result states that the fundamental tone of a free plate under tension is 
sharply bounded above by the fundamental tone of the ball of same volume, which 
is analogous to the Szgo- Weinberger result for free membranes. 

The circular plate is also significant because it is a rare example where we are 
able to find the form of its eigenf unctions. The eigenvalue equation for the free 
plate is not separable in rectangular coordinates due to the cross-term Uxxyy in 
AAm, and so we cannot easily find the eigenfunctions for the rectangular plate. 

Problems with the bi-Laplacian tend to be more difficult than their second-order 
counterparts, because the theory of the bi-Laplacian is not nearly so well developed. 
For example, the maximum principle fails for the bi-Laplacian, and solvability of 
the biharmonic equation in Lipschitz domains with Neumann boundary conditions 
was established only a few years ago by Verchota [21 . 

Other notable mathematical work on plates includes Kawohl, Levine, and Velte 
[8] , who investigated the sums of low eigenvalues for the clamped plate under tension 
and compression, and Payne [TBj , who considered both vibrating and buckling free 
and clamped plates and established inequalities bounding plate eigenvalues by their 
(free or fixed) membrane counterparts. 



2. Formulating the problem 

In this section we present the mathematical formulation of the free plate problem, 
summarizing known facts about existence of the spectrum, boundary conditions, 
and regularity of solutions as proved in [5]. We consider dimensions d > 2; the 
elementary one-dimensional case (free rod) is discussed in [6j Chapter 7]. In this 
paper we consider only regions $7 = IB(^), balls of radius R\ we write B = B(l). 

The eigenvalues of the free plate under tension with shape B(i?) are characterized 
by the Rayleigh quotient 

Here \D'^u\ = {Y.jk^l,xJ^'^ is the Hilbert-Schmidt norm of the Hessian matrix 
D^u, and Du denotes the gradient vector. The parameter r > has the physical 
meaning of the ratio of tension to fiexural rigidity. Note in particular that we have 
the Hessian term |Z3^up, involving all second derivatives, rather than only those of 
the Laplacian term (Aw)^; see [321 P- 228]. 

From this quotient, we obtain the Euler-Lagrange partial differential equation 

(2) AAu - rAu = wu, 

where uj is the eigenvalue, together with the natural (i.e., unconstrained or "free") 
boundary conditions on dM{R) [3 Proposition 7]: 

(3) Mu u,.,, = at r = i?, 

(4) Vu TUr - ^As (ur - ^) - i^u)r = at r = R. 
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Here Ag is the angular part of the Laplacian, which in two dimensions is simply 

d'^/de'^. 

The spectrum consists of nonnegative eigenvalues of finite multiplicity 

Q — UJq < UJi < 102 < ■ ■ ■ ^ GO. 

The corresponding eigenfunctions are smooth on B(i?). The eigenfunction corre- 
sponding to Wo = is the constant function. 

We call LOi the fundamental tone and a corresponding eigenfunction ui a funda- 
mental mode. 

The radial part of the eigenfunctions for the ball will be obtained in terms of 
ultraspherical Bessel functions. In order to state the main results of this paper, we 
will need to define these functions. The ultraspherical Bessel functions ji{z) of the 
first kind are defined in terms of the Bessel functions of the first kind, Jv{z), as 
follows: 

ji{z) = z^^Js+iiz) 

where the parameter 

d-2 



depends on the dimension d. This notation suppresses the dependence of the ji 
functions on the dimension d, which causes no problems since the dimension d> 2 
is fixed. Ultraspherical modified Bessel functions ii{z) of the first kind are defined 
analogously: 

ii{z) ^ z^-''Is+i{z), 
where Ii, is the modified Bessel function of the first kind. 

3. Main Results 

The aim of this paper is to find the eigenfunctions of the free plate problem 
on the ball and identify the fundamental mode. We will only consider r > 0, 
corresponding to the free plate under tension. 

The first result identifies the form of the eigenfunctions of the ball of radius 
R>0. 

Proposition 1. (Eigenfunctions in spherical coordinates) Let r > and uj be any 

positive eigenvalue of the free ball'B(R); that is, uj is an eigenvalue of AAu — tAu = 
LOU under boundary conditions ^ and Then the corresponding eigenfunctions 
can he written in the form Ri{r)Yi{9), where Yi is a spherical harmonic of some 
integer order I and Ri is a linear combination of ultraspherical Bessel and modified 
Bessel functions, 

Ri{r)= ]i{ar/R)+^ii{hr/R). 

Here the positive numbers a and h depend on r and lo by h"^ — ~ R^t and 
a^b^ ~ R'^LO, and ^ is a real constant given by 

^ bH'{{b) ■ 

This result gives us the forms of eigenfunctions given an eigenvalue w. However, 
we do not know the eigenvalues, and instead wish to find them. The next result, a 
consequence of the proof of Proposition [H allows us to compute the eigenvalues of 
the ball: 
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Proposition 2. (Eigenvalues) Fix r > 0. Then a real number uj is an eigenvalue 
of the free plate of shape B(i?) if and only if uj = or there exists an integer I > 
and positive real constants a, b such that a^b^ = R*uj, b'^ — ~ R^t, and 

Wi{a) :=aV(a) (^-a^bi'i{b) + l{l + d - 2)(^bi'i{b) - ii(b)'^^ 

- b\'l{b) (ab^j'iia) +l{l + d- 2)(a3l{a) ~ j,(a))) = 0. 

Remark. This proposition gives us a way to calculate eigenvalues numerically. Given 
T and a nonnegative integer I, the roots of Wi{a) will determine eigenvalues by the 
relation uj — a^la^ + t). 

We are particularly interested in identifying the fundamental tone and mode of 
the ball, because the proof of the free plate isoperimctric inequality in [5 uses trial 
functions constructed from that fundamental mode. Thus the following theorem is 
essential to the proof in [5| . 

Theorem 3. (Fundamental mode of the ball) For t > 0, the fundamental modes 
of the ball B(i?) can be written as linear combinations of 

Mr J) = (ii{ar/R)+^ii{br/R))Yi{e), 

with a, b, 7 real constants, with a and b positive and depending on t and iO\ as 
follows: b^ — = R^t and a^b^ = R'^oJi, and 7 given by 

^ bH['{b) ■ 

In particular, in dimension 2, the fundamental modes are linear combinations of 
ui = ( Ji (ar/i?)+ 7/1 (6r/i?)) 

Remark. When r = 0, linear functions are eigenfunctions with eigenvalue zero, and 
so the zero eigenvalue is d + 1-fold degenerate; the fundamental modes will then 
involve higher-order spherical harmonics. 

This model considered with negative values of the tension parameter r corre- 
sponds to the free plate under compression. In this case, the Rayleigh quotient Q 
yields both positive eigenvalues (corresponding to vibrational modes) and negative 
eigenvalues (corresponding to buckling modes) . The forms of the eigenfunctions in 
the buckling modes differ from those of Proposition [TJ Furthermore, for t < the 
fundamental modes can involve higher-order spherical harmonics. 

4. Ultraspherical Bessel Functions of the First Kind 

We must examine properties of d-dimensional ultraspherical Bessel functions, for 
they provide the radial part of the eigenfunctions of the ball. For more information 
on Bessel functions, see p. 358-389]. For more information on spherical and 
ultraspherical Bessel functions, see [U p. 437-455] (d = 3 only) and [12] (all d>2). 

The Bessel ji (z) functions defined previously solve the ultraspherical Bessel equa- 
tion, 

(5) zV + {d- l)zw' + (^2^ -l{l + d-2)^w = 0. 
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The modified Bessel ji{z) functions solve the modified ultraspherical Bessel equa- 
tion, 

(6) z^w" + {d- l)zw' - (^z^ + l{l + d-2ijw = 0. 

Recurrence Relations and power series. The Bessel functions J^, and Ii, have 
a number of useful recurrence relations; see, for example, jTJ p. 361, 376]. The 
ultraspherical Bessel functions have similar recurrence relations, all of which follow 
from the definition and application of the corresponding ordinary Bessel recurrence 
relations. 

(7) - — ^ ji-ii^) + ji+ii^) 

(8) jliz) = ljiiz)-Ji+i{z) 

(9) =,n^,{z)-^-±^,n{z) 

(10) ii(z) = ii_i[z) - ii+i[z} 

z 

(11) i'liz) = -ii{z) + ii+iiz) 

z 

Note that if we take d — 2, each of these simplifies to the corresponding relation 
for Bessel functions. 

We also have recurrence relations for the second derivatives: 

(12) j'/{z) = - l) Jiiz) + ^Ji+i{z) 

(13) ,;'(;,)= + ,,(,)_ ^^,^^(^). 

Again, when d — 2 each recurrence relation simplifies to its two-dimensional analog. 

We may also write a power series for the ultraspherical Bessel functions ji{z) 
and ii{z) using the series for the corresponding Js+i and Is+i- 



(14) Ji{z) = Yl 



(15) ^i{z)^J2 



,4fc!r(fc + f + 



By examining the power series (|15p . it is immediate that ii{z) and its derivatives 
are all positive on (0, oo). Since the terms of the power series for ji and ii are the 
same up to a sign, we also have that the derivatives of ji are dominated by those 
of if. 



(16) jj^H^) <»P(^) form>0, z>0, 

with equality only at z = 0. 

Other needed facts. To identify the fundamental tone of the circular plate, and 
to prove the isoperimetric inequality for the free plate in [5], we will need several 
facts about Bessel functions and their derivatives. We begin with a result on the 
zeros of the j'i{z). 
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Proposition 4 (L. Lorch and P. Szego, [T2])- Let pi denote the kth positive zero 
of jl{z). Then for d > 3 and I > 1, 

l{d + 2l){d + 2l + 2) ,2^,,,^„,, 
j:^^^ <ip,,) <l{d + 2l). 

In particular, for pi^i the first zero of j[, we deduce 

d < pI -^ <d + 2. 

This inequality holds for all d > 2. 

Recall s — (d — 2)/2. This first lemma below gives us information on the signs 
of the Bessel ji functions, and will be used in to help prove the free plate 
isoperimetric inequality. 

Lemma 5. The functions ji and Jg+i have the same sign. In particular, for I = 
1, . . . ,5 and any d>2, we have ji{z) > for z < pi^i- 

Proof. The first statement is immediate from the definition of the ultraspherical 
Bessel functions. For the second statement, we appeal to established facts of Bessel 
functions. We write for the first nontrivial zero of the Bessel function Ji{z). It 
is well known that Ji{z) is positive on {0,ji.i) and the zeros ji^i are increasing in 
I for I > 1. Because Ji(z) = at z = and ji.i with no zeros between, we have 
the same for ji{z) and thus the first root of pi.i, lies between and ji^i. 

Therefore for any d > 2 and any I > 1, we have Ji{z) > and hence ji{z) > on 
(0,Pi,i]- □ 

Lemma 6. We have j[ > on (0,pi^i). 

Proof. This follows from the observation that ji{z) > on (0, ji^i) with ji{0) — 
and the definition of pi i as the first zero of j'i(z). □ 

Lemma 7. We have ji^ > on (0,pi^i]. 

Proof. Let p2,i denote the first zero of By Proposition 01 pi^ < d + 2 and 

2 ^ 2(rf + 4)(d + 6) 

^^2.1 > ■ 

Thenpli ~pIj > + 8d + 28)7(10 + d) > 0, so that > on (0,pi,i]. □ 
Lemma 8. We have j'( < on (0,pi,i]. 
Proof. We see that 

f;{z) = ^j2{z)-n{z) hym 

= -\h{z) - Uz) by dSl) with / = 2. 

Since both j2 and ^2 are positive on (0,pi,i] by the previous lemmas, we obtain 
j'l > on that same interval. □ 

The next two lemmas are needed for the proof of the free plate isoperimetric 
inequality in p; but are not used in this paper. We include them here since both 
follow from properties of Bessel functions established in this paper. 

Lemma 9. We have j'f^ > on (0,pi^i]. 
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Proof. We have by ^ that 
and so 

• / ^ 2(d- 1) . , , d2 - 1 , ^ 
Z Z'^ 



by (ini) with I = 1 and / = 2, and (jH]) with I = 2. When d = 2, this becomes 
(18) jf)(z)=('l-A')ji(z)+('i^--')j2(z) 



z 



by with ; = 2. For any d, (d?]) gives us 

(19) jf{z)^^—^3-,{z)+('^^^-\\3^{z) by with / = 2 



(20) ^ ( i^_i)j3(z) + i-ij4(z) byOwith; = 3 



(21) = ( i^^^ - J4(z) + fl - J5(-) by with / = 4. 



z^ 



When c? = 2, then the first term of (fTO|) is nonnegative on (0,^1^] by Lemma [T] 
The function jj, is positive on (0,pi.i] by Lemma[5J note that since c? = 2, we have 
j3(z) = J'i{z). Thus we have i'^^ [z) > when z G (0, -y/S] n (0,^14]. However, 
Pi,i > V^, so we have only estabhshed positivity on (0, \/3]. 

To estabhsh positivity on (-\/3,Pi,i] we turn to (fT8|) . The first term is certainly 
positive on {y/3,pi,i]. The second term is positive when both J2 > and z < \/6. 
Because ~ 1.84 for d = 2, we have pi^i < and we are done. 

When d — 3 and d = 4, we again examine (IT^ . Then Lemma [S] together with 

the argument above give us > on (0, Vd^ — 1) n (O,^!^]. By Proposition |4] 
we have < \/d + 2, which for d = 3 and 4 is less than \/ (P = 1 , thus proving 
the lemma for these d. 

For dimensions d > 5, we turn to pop . The second term is positive on (0,pi.i] 
for all d > 4 by Lemma [5) Since pi_i < \/(ri~2 and \^d + 2 < vTs for d < 13, we 
conclude ji^\z) > on (0,pi,i] for d = 5, . . . , 13. 

Finally, suppose d > 14 and z G (0,pi.i]. If z e (0, \/T5], then /i^^z) > as 
above. If z > \/T5, then we examine ([2T|) . Here the first term is nonnegative on 

[■\/T5,pi_i]. The non-Bessel factor of the second term is positive on (0, \J \{d + 6)] 
and hence on (0,^1^]. □ 

The next lemma provides some bounds on ultraspherical Bessel functions that 
will be needed in |5] and will not be used in this paper. 

Let dfc denote the coefficients of the series expansion for j"(z), so that 



00 00 
Y^i-lfduz^^-' and i'i{z) = Y,dt.z^^- 

k=l k=l 
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by (HH) and ((T5l), where 

dk 



2k + I 



(fc-i)!r(fc + i + d/2) 

Lemma 10. We have the following hounds: 



2l-2fc-d/2 



-diz + ^22^ > ii{z) 
6 

diz + -d2Z^ > i'l{z) 
5 



for all z G 



0, V3(d + 2)/(d + 5) 



/or aZ/ z G 



0,\/3 



Proof. Let 



2/S + 3 



Cfc 



dfc ^ 2fc(2fc + l)(2fc + rf + 2)' 



It is easy to show that Ck is decreasing for A: > 1. 

We use the series expansion to first prove the following upper bound on 
for z > 0: 



{-d,z + d^z^) - f{(z) = 

A.— 3 
oo 

= ^ (1 - Ckz'^)dkZ 



k=3 
fcodd 



> (1 - ciz2) J2 d, 



kZ 



^2k~l 



fc=3 
/codd 



since Cfc is decreasing in k. Hence {~diz + d2Z^) — ji{z) > for all z with < 
z < l/VcT = ^/6id + A)/5, which is a larger range even than claimed in the first 
estimate in the lemma. 

For we take a slightly different approach. We will show that on [0, \/3], 



^ oo 
-d2Z^ > ^42' 



fc=3 



and thus 
(22) 

On [0, -\/3], note that 

00 00 



diz + -d2Z^ > ii{z). 
5 



2fc+ 1 



k=3 



k=Z 



(fc-i)!r(fc + i + d/2) 



2-d/2 



2fc-l 



2fc + 1 



2fc-4 



fe=3 

since z < Vs 



{k - iy.{k + d/2)r{k + d/2) 
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^ 2-rf/2 rz\3^ 2fc + l /3y~' 
-r(3 + d/2)V2/ ^ (fc- l)!(A: + d/2) U/ 
since r(2:) > 1 and is increasing on [2, oo), 

- 5 ^ fc! U; 

fc=2 ^ ^ 

by the definition of ^2 and taking fc i— > fc + 1 
by tfie power series for 

5 

Thus we have obtained our desired bound on i'/. □ 



5. Ultraspherical Bessel functions of the second kind 

Each of the Bessel equations ([5|) and ([1]) is a second-order differential equation, 
and so has another set of solutions — namely, Bessel and modified Bessel functions 
of the second kind. However, these functions are singular at the origin. By the 
regularity of plate eigenfunctions, either these singular solutions do not appear in 
the eigenfunctions, or they appear in a linear combination such that the singular 
terms cancel. Lemma which appears below, states that in fact there is no 
nontrivial linear combination that meets the smoothness condition. 

Ultraspherical Bessel functions of the second kind solve ([5]) and are defined by 

miz) = z-'Ns+iiz) 

with s — {d — 2)/2. Here N^{z) denotes a Bessel function of the second kind of 
order v. Each Ny{z) is linearly independent of Jv{z) (see, for example, [U p. 358]), 
so ni{z) is linearly independent oi ji[z). The functions Ni{z) are often written as 
Yi{z); we use Ni to avoid confusion with the spherical harmonics Yi{9). 

Ultraspherical modified Bessel functions of the second kind solve (O and are 
defined by 

ki(z) = z-'Ks+i{z) 
d-2 

with s = , 

2 ' 

where K^{z) denotes a modified Bessel function of the second kind of order v. The 
ki{z) are linearly independent of the ii{z). 

Lemma 11. Let a, b be positive constants, a < b. 

For d > 2 and all integers I > 2, there is no nontrivial linear combination 

R{z) = Ani{az) + Bki{bz) 

so that R{z) is smooth at z = 0. 

For d > 2 and I — 1, there is no nontrivial linear combination 

R{z) = Ani{az) + Bki{bz) 

so that R{z) is smooth at z — with i?(0) — 0. 
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For d> 2 and I = 0, there is no nontrivial linear combination 

R{z) = AnQ{az) + Bko{bz) 
so that R{z) is smooth at z — with R'(0) — 0. 

The proof of the lemma is technical and so has been omitted; it can be found in 
full in % Chapter 4]. 

6. Proof of Proposition [T] 

The full set of eigenfunctions for the circular free plate under tension will be found 
exactly in terms of Bessel and modified Bessel functions and spherical harmonics. 

We will focus on the unit ball, since the solution of our eigenvalue problem for 
any ball B(i?) can then be obtained by scaling, as follows. If u{x) is an eigenfunction 
of the unit ball with eigenvalue uj and tension t, then u{x) — u{x/R) is an eigen- 
function of the ball B(i?) with eigenvalue i?^a; and tension t/R?; see [F, Lemma 
24] . We will establish Propositions [T] and [21 which give us the eigenfunctions and 
eigenvalues of the ball. Recall that Proposition [T] states that all eigenfunctions will 
be of the form Ri{r)Yi{9), where Ri is a linear combination (depending on r) of 
ultraspherical Bessel and modified Bessel functions of order I, and Yi is a spherical 
harmonic. 

The eigenfunctions of the free disk with zero tension are treated in [TTl Chap- 
ter 2], along with treatment of the free disk with clamped and simply supported 
boundary conditions. In addition to finding the eigenfunctions, Leissa provides 
some numerical computations of eigenvalues and presents experimental data on 
such plates. 

Spherical harmonics. In the case where is the ball, it is natural to consider 
spherical coordinates. Let r be the radius and 9 be the remaining angular informa- 
tion. Consider Laplace's equation A/ = 0, with / a function on M.'^. The Laplacian 
can be written in spherical coordinates as 

A = drr + ^y^9r + ^Ag, 

where we give the name A5 to the angular part of the Laplacian. Separating 
variables so that / = R{r)Y{9), we obtain 

R" + ^^R' -litt^^R = and AsY = -l{l + d - 2)Y. 

The parameter I appearing in the separation constant 1(1 + d — 2) must be an 
nonnegative integer in order for solutions to exist. The solutions to 
AsY = —l(l + d — 2)Y are called the spherical harmonics. For each I, we choose 
a spanning set {Yi} of such solutions that are orthonormal in the L^((9B) inner 
product. Because the eigenvalues are real, the Yi may be chosen to be real-valued. 
However, they are traditionally chosen to be complex- valued, and so will be treated 
as possibly such in the proof of Theorem |3l 

Factoring the eigenfunction equation. The proofs of Propositions [T] and [2] 
involve a factoring of the eigenvalue equation akin to factoring of the characteristic 
equation of an ordinary differential equation. It is notable that while the eigenvalue 
equation is not separable in rectangular coordinates, the factoring together with the 
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commutativity of the spherical part of the Laplacian As with the Laplacian A gives 
us a sort of separability in spherical coordinates. 

Proof of Proposition QJ We first show that eigenfunctions can be written as a prod- 
uct of a radial function with a spherical harmonic, and then give the exact form of 
the radial part. 

Write A := — tA. From Section 2, we know each eigenvalue w has finite mul- 
tiplicity, and so the corresponding space of eigenfunctions X^^ is finite-dimensional. 
Because As is independent of r, it commutes with the Laplacian A and hence with 
our operator A = A^ — rA. Thus As maps X^j into itself. The operator As is sym- 
metric, and so is diagonalizable on the finite-dimensional space X^j- Thus we have 
that A and As are simultaneously diagonalizable. The eigenfunctions of As on dM 
are the spherical harmonics; on B the eigenfunctions have the form R{r)Yi{6). Wc 
can therefore choose our eigenfunctions of A to have this form. 

To find the precise form of R, we factor the eigenvalue equation ([2|), obtaining 

(23) (A + a2)(A-62)u = 0, 

where a and b are positive real numbers satisfying b"^ = a? -\- t and lu — a?' {a? + r) . 
That is, = v^(t/2)2 +a;-T/2 and 6^ = ^(t/2)2 +u + t/2. The eigenfunctions 
u will then be linear combinations of the solutions vi and V2 of each factor: 

(24) (A-|-a^)wi=0 and {A-b^)v2^0. 

Each of these is separable in spherical coordinates, with angular equation 
AsY — + d — 2)Y for some nonnegative integer The radial equation for vi 
is a rescaling of the ultraspherical Bessel equation ([5]) with order / and the radial 
equation for V2 is a rescaling of the ultraspherical modified Bessel equation ^ with 
order Z, hence 

vi = (^Aji, (ar) + Bni, (ar)) Yi, and V2 = (ci;, (br) + Dki, (6r)) Yi, , 

for some nonnegative integers Zi, I2 and real constants A, B, C, and D. From the 
diagonalization argument above, we know u = R{r)Yi{6), so all the orders must 
agree: I = li = h- Thus solutions of the eigenvalue equation © have the form 

uix) = Rir)T{§) - (^Ajiiar) + Bniiar) + Ciiibr) + Dki,ibr)^Yii9). 

However, the eigenfunctions are smooth on B. The spherical harmonics Yi have no 
radial dependence; thus we must have the radial part R{r) be a smooth function of 
r G [0,00). When I = 0, the spherical harmonic Yq is constant, and we must also 
require R'{0) = in order for u to be smooth. When 1 = 1, the spherical harmonics 
Yi can be given by Xi/r, where Xi are the coordinate functions. Then along the 
Xi-axis, R{r)Yi{0) = R{r)xi/r = R{r) sign{xi). This function is continuous only if 
R{r) vanishes at r = 0. By Lemma [Tl] there is no nontrivial linear combination of 
Bessel functions of the second kind which satisfies these conditions; thus B and D 
are both zero. Denote C/A by the constant 7; then we have 



uix) = R{r)T{9) = (^ji{ar)+jii{br)yii§). 
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The constant 7 must be chosen so that u satisfies the natural boundary condition 
Urr = at r = 1; hence we have 

and so 7 is real-valued. □ 

7. Proof of Proposition [2] 

We prove Proposition [5] for the unit ball; the general result follows from scaling 
as noted previously. Most of the work has already been done for us in the proof of 
Proposition [TJ 

Proof of Proposition [H It is immediate from the Rayleigh quotient Q that uj — 
is an eigenvalue with corresponding eigenfunction u — const. 

Fix r > 0. Now suppose w ^ is an eigenvalue. Then by Proposition [T] we must 
have an associated eigenfunction u of the form (^ji{ar) + 'yii(br)j Yi for some index 

I and positive real numbers a and b with b = y/a'^ + t and uj — a^b^ . 

All eigenfunctions u satisfy the natural boundary conditions, Mu = and 
Vu = on dil, as given in (|3]) and Q for the ball. Since by Proposition [T] all 
eigenfunctions are linear combinations of ji{ar)Yi{9) and ii{br)Yi(9), we must have 
some nontrivial linear combination satisfying the homogeneous linear equations 

jMu = 
[Vu = 0. 

Thus we must have that the determinant 

^ ^ \Vji{ar) Vii{br))^^^ 

vanishes. From the first natural boundary conditions for the ball given in ([3]), we 
have 

(25) Mji{ar) = a^f/iar) and Mii{ar) = b^i'l{br). 

The ji{ar) and ii{br) are rescaled ultraspherical Bessel and modified Bessel func- 
tions, so by the factorization ([24|) . we have 

Aji{ar)Yi{9) ^ -a^ji{ar)Yi{9) and Aii{br)Yi{9) ^ b^ii{br)Yi{9). 

Then noting r = 1 on dM, the "V boundary condition terms from in ^ can be 
rewritten as follows: 

Vjiia) = raj'iia) + l{l + d - 2) (a][{a) ~ j,(a)) + a^j',{a) 
Vii{b) = Tbi'i{b) + l{l + d - 2)(bi[{b) - ii{b)^ -b^i'i{b) 

Combining the above with ((25l) and substituting t = b^ — a? , we find 
Wi{a) ^ a2j;'(a)( - a'^bi[{b) + l{l + d- 2){bi'i{b) - i,(5))) 

(26) - b^'ib) [ab'j'iia) + lil + d- 2)(aj,'(a) ~ Ma))) . 
Because to is an eigenvalue, the determinant Wi{a) vanishes, as desired. 
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Now suppose we have a number w with some integer / > and positive real 
numbers a, b, satisfying u — a^b^, b^ — = t, and Wi{a) = for some nonnegative 
integer I. Then the function u = (ji{a)+jii{b)^Yi satisfies the boundary conditions 

Mu = and Vu = for r = 1. Furthermore, u satisfies the eigenvalue equation 
(121) with eigenvalue a^fe^ = w, and so u) is an eigenvalue. □ 

8. Proof of Theorem [3] 

In this section, we identify the fundamental mode of the ball for positive tension, 
proving Theorem |3l 

The proof will have two parts. First we show that for any radial function R{r), 
the Rayleigh quotient Q[RYi] is minimized when / = 1, among alH > 1. Then we 
show that of all nonconstant eigenstates with I — and ^ = 1, the lowest eigenvalue 
corresponds to I = 1. Note that when / = 0, the spherical harmonic Yq is the 
constant function, and so I = corresponds to purely radial modes. 

Proof of Theorem 3. [Part 1.] We will show that for any fixed smooth radial func- 
tion R, the Rayleigh quotient Q[RYi\ is an increasing function in I for all Z > 1. 
Then by the variational characterization of eigenvalues, we see that the lowest 
eigenvalue corresponding to an eigenfunction with angular dependence (i.e., / > 1) 
occurs when 1 = 1. 

Considering the numerator and denominator separately, we will use the (9IB)- 
orthonormality of the spherical harmonics to simplify the angular parts of the in- 
tegrals. 

The denominator of our Rayleigh quotient is, for u = RYi^ 

[ \RYi\^dx= [ R^r'^-^dr, 
Jb Jq 

and so is independent of L So it suffices to show that the numerator is an increasing 
function of I for I > 1. 

Recall the numerator of the Rayleigh Quotient is 

N[u] = I iD^ul"^ +t\Du\^ dx. 

We use the following pointwise identity to rewrite the Hessian term: 

(27) \D^u\^ = i (^A{\Du\^) - D{Au) ■ Du - D{Au) ■ Du^ 

Because our region D, is the unit ball, we may use spherical coordinates, noting 
^ = Ur- We then write the gradient as 

(28) Du = Urf + -VsM, 

r 

where f = x/r \s the unit normal. Note that ^Vs is the surface gradient gradgg of 
the ball, and ^As is the Laplacian on the boundary of the ball. We write div^a 
for the surface divergence of the ball; then ^Asf = div^Bgradg^/. 
Then by the Divergence Theorem on 9B, we have for any function /, 

(29) / ^AsfdS^ f divaB(gradaB/)d5-0. 
Jam''' Jam ^ ' 
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Exploiting orthonormality of the Yi, we see that 
(30) / \\7sYi\^dS^l{l + d-2)^:k. 



Thus when u = RYi, we use (|27p to rewrite the integral of the Hessian term as 
follows: 

i J (^A{\Du\^) - D{Au) ■ Du - D{Au) ■ Du^ dx 

r^- r^f d-1 k \ ^ ^(_ d-\_ k_\ , 
— Du ■ D I Urr H Ur -U | — Du ■ D | Urr H Ur -U \ dx 



r r- 
(noting AsYi = -kY) 



d-l k \ ^ d-l_ k_. , 

Du ■ D Urr H Ur rrM — Du ■ D Urr H Ur ttU dx , 



with this last by noting that (|29l) gives us 

/ ^As (iRr\Yi\' + ^\ysYiA dx = 0. 
Expanding the integrands, our integral of the Hessian term simplifies to 



Then integrating the above over IB using (pOl) and the orthonormality of the Yi , we 
obtain 



/ \D^ufdx 



V, /,n9 2fc + d-l, 6fc , k(k-d + A) A . . , 
I {R"Y + [R'f - —RR' + '-R^ r'^-^ dr 







V, d-l, 2k/ , 3 \ k(k-d- 1/2) ,\ ^ , 
(R'T + —riR'f + -j[rR'~-R] + ^ -^^R' r'^' dr 



with this last equality by completing the square. 

We now examine the gradient term in the numerator of the Rayleigh quotient: 



/ \Du\^dx= ( + 4|Vswn dx 



(RTlYif + ^l^sYif^ dx 



again by ([50]) . 
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Combining these results, the numerator of the Rayleigh quotient can be now 
written with all fc- dependence (and hence Z-dependence) explicit: 

N\u]= I [^[rR'--R\ +'-^ R:^ + t'-^ | r^-'dr 




Recall r > 0; then the above is increasing with k ioi k > d + 1/2. Recall 
k = l{l + d — 2) is increasing as a function of I; then all terms involving k are 
increasing functions of I for I > 2. If Z = 1, the expression k{k — d — 1/2) becomes 
— 3(d— l)/2, which is negative for all dimensions under consideration. If / = 2, we 
find k{k — d — 1/2) = 2d{d — 1/2) > 0. Thus each term involving I is increasing as 
a function of I for alH > 1. 

Thus for any fixed radial function R, the numerator N[RYi] is an increasing 
function of / for Z > 1 . 

[Part 2.] We now show that the lowest eigenvalue corresponding to an eigen- 
function of the form ui = (ji{ar) + 'yii{br)^Yi with Z = 1 is less than the lowest 

positive eigenvalue for / = 0. 

Let denote the first positive zero of j[ (a) . Recall from Proposition [2] that uj 
is an eigenvalue if and only if we have some integer I > and can write u = a'^b'^ 
for positive constants a and b such that b^ — a? — t and Wi{a) = 0. The parameter 
T is positive, so w = a^(a^ +t) increases with a. Therefore, to show that the lowest 
nonzero eigenvalue corresponds to Z = 1 and not Z = 0, we show that the first 
nonzero root of Wi (a) is less than the first nonzero root of Wq (a) . 

First we consider I — 0. Here A: = 0, so we look for solutions to: 

= Wo (a) 

= a^f^{a)[- a%r',{b)) b\'^{b)(ab^ Ua)) 
= a%j[{a)ii{b) + ab'^i'^{b)ji{a), 

by ([8]) and (fTTj) . The functions jo(6) and «o(6) are positive for & > 0, as noted earlier 
in this paper, by the power series expansion (flS)) . Similarly, ji(a) and j'i{a) are 
positive on (0, Ooo) by Lemma [SI and so Wo(a) > on (0, floo)- 
Now consider I = \. The constant fc = d — 1, so we have 

Wi{a) - a2jY(a)( - aHi[{b) + {d-l){bi[{b) - zi(6))) 

(31) - bh'l{b)(ab'j[{a) + {d-l){a3[{a)-j,{a))) 
= a'b3';{a)(^- a\\{b) + {d- l)^2(fe)) 

(32) -abX{b){f3>)-{d-l)n{a)), 

by the Bessel identities ([8|) and (fTTj) . As a 0, the first term in ([32]) behaves like 

o^VTji{a)(^[d — 1)12(^/1-)^ 

and so its sign is determined by that of j" (a)- By Lemma HI i"(a) is negative for 
all a e (0, Ooo]- Hence the first term of Wi{a) is negative near a — Q. 
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Similarly, as a — >■ 0, we see ^2(0) — > 0, and so the second term behaves like 

which is negative by Lemma |51 Therefore Wi{a) < near = 0. 

At a — Coo, we have Ji(aoo) = by definition of aoo! also note ji(aoo) > and 
ji(aoo) < by Lemma [5] and ISl Write 600 — vaoo+^- By Proposition |31 we 
have > d > d - 1 for d > 3; for d = 2, we have 1.84^ > d - 1. Thus 

^00 — {d — 1) > 0, and so (PT|) gives 

W^i(aoo) = -a^iJ'(aoo)(6oo(aL - {d - l))i[{boo) + {d-l)ii{b^)) 

+ bli';ib^)id-l)nia^). 

Hence both terms in Wi(aoo) and since Wi is continuous, it must have a zero in 
(OjOtx))- Thus the lowest nonzero root of Wi{a) occurs when / = 1, not I = 0, and 
so the lowest eigenvalue oj — a{a'^ + t) occurs when 1 = 1. □ 
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